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ABSTRACT 

We have obtained X-ray observations of the bipolar planetary nebulae (PNe) NGC 2346 and 
NGC 7026 with XMM-Newton. These observations detected diffuse X-ray emission from NGC 7026 
but not from NGC 2346. The X-ray emission from NGC 7026 appears to be confined within the bipolar 
lobes of the PN and has spectral properties suggesting a thermal plasma emitting at a temperature of 
l.ltoi X K. The X-ray spectrum of NGC 7026 is modeled using nebular and stellar abundances 
to assess whether a significant amount of nebular material has been mixed into the shocked-wind, 
but the results of this comparison are not conclusive owing to the small number of counts detected. 
Observations of bipolar PNe indicate that diffuse X-ray emission is much less likely detected in open- 
lobed nebulae than closed-lobed nebulae, possibly because open-lobed nebulae do not have strong fast 
winds or are unable to retain hot gas. 

Subject headings: planetary nebulae: general — planetary nebulae: individual (NGC 2346, NGC 7026) 
— X-rays: ISM — stars: winds 



1. INTRODUCTION 

Planetary nebulae (PNe) consist of the stellar mass 
lost by their low- or intermediate-mass stellar progen- 
itors during late evolutionary stages. As a star evolves 
past the asymptotic giant branch (AGB) phase, the mass 
loss process changes from a copious slow AGB wind to 
a tenuous fast wind and the interaction between the 
fast wind and the AGB wind forms a PN (yKwok 1983.1 . 
iFrank fc Mellemal l)1994j) have demonstrated that this 
wind-wind interaction can produce the commonly ob- 
served elliptical and bipolar morphologies of PNe, al- 
though the bipolar morphology requires that either the 
AGB wind or the fast stellar wind is anisotropic and has 
a polar angle dependence. 

The geometry of the fast wind has often been assumed 
to be isotropic, but the presence of bipolar jets and col- 
limated outflows in some PNe indicates that the fast 
wind may span a range of geometries. For wind ve- 
locities greater than ~300 km s~^, the interaction of 
the fast stellar wind with the AGB wind will shock-heat 
the wind material to temperatures greater than 10® K. 
This shocked fast stellar wind is too tenuous to pro- 
duce appreciable X-ray emission, but at its interface with 
the dense nebular material, heat conduction lowers the 
temperature and mass evaporation raises the density of 
the shocked fast wind, producing optimal conditions for 
X-r ay emission (|Weaver et aLllTOT^ iZhekov fc Perinottol 
Il996il Since the hot gas in a PN is associated with 
shocked fast wind, the distribution of X-ray emission 
may be used to infer a directional dependence for the 
fast wind. 

Chandra and XMM-Newton observations of ellipti- 
cal PNe have revealed diffuse X-ray emission confined 
within their innermost nebular shells , e.g., NG C 2392 
lIGuerrero et aIJl2005all and NGC 6543 ijchn et aIJl200ir . 
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The plasma temperatures implied by their X-ray spectra 
are a few times 10® K. In the case of NGC 6543, the X-ray 
emission is well resolved, and a clear limb-brightening is 
observed. Of the ten elliptical PNe observed by either 
XMM-Newton or Chandra, six show detectable diffuse 
X-ray emission l|Guerrero et alJl2005b|) . 

Chandra has observed many bipolar PNe. While two 
of them, Mz 3 and NGC 7027, exhibit diffu se X-ray emis- 
sion confined within the bipolar lobes feastner et alJ 
l200li r2003). the other bipolar PNe (e.g., MyCnlS, M 1- 
16, Ml-92, M2-9, OH 231.8+4.2) were not detected. 
Most of these non-detections are understandable because 
they are associated with relatively young/exotic bipo- 
lar PNe that either have large circumstellar absorption 
or have coUimated outflows with low wind/outflow ve- 
locities, <300 km s^^. In order to better characterize 
diffuse X-ray emission from bipolar PNe, we have used 
XMM-Newton to observe two carefully chosen bipolar 
PNe with contrasting lobe morphologies, NGC 2346 and 
NGC 7026. 

NGC 2346 has a butterfly morphology with two open 
bipolar lobes and a bright equatorial waist that co ntains 
large amounts of molecular gas and dust (Walsh ~t alJ 
IT99ltlZnckerman fc Gatle vl988: 's7ei~aLl 20041. Rs cen- 
tral s tar is a white dwarf with an A5 companion ( AIcnde3 
^2l)j forming a si ngle-line spectroscopic bin ary with a 
period of ~16 days l)Mendez fc Niemelalll98lD . Its nebu- 
lar gas shows a bip olar velocity field with a dynamic ag e 
of 3,500-4,700 yrs ('Wa,l sh et a].HT99TtlAria,s et a,].ll2001^ . 
NGC 2346 can be described as an evolved bipolar PN. 

NGC 7026 has closed bip olar lobes emanating from 
a bright ring-like waist l|Solf fc Weinberger! Il984() . 
and the lobes appear to ha ve m ultiple compo nents 
(jCuesta. Phillips, fc Mampasol Il996t iLoped l2003|). Its 
central star is a hy drogen-deficient WC star ijAlleJ 
119761 lAcker fc Nein er 2003). The dynamic age derived 
from the nebular expansion, < 1,000 yrs, suggests that 
NGC 7026 is a young bipolar PN (ISolf fc Weinberger! 
ri98l . 

The XMM-Newton observations of these two PNe have 
detected diffuse X-ray emission from NGC 7026, but not 



2 



Gruendl et al. 



from NGC 2346. In this paper, we describe these X-ray 
observations and their reduction in §2, report our anal- 
ysis of the observations in §3, and discuss their imphca- 
tions in §4. 

2. OBSERVATIONS 

NGC 2346 and NGC 7026 were observed with XMM- 
Newton in Revolution 876 (Obs ID 0200240501) on 2004 
September 20 and Revolution 825 (Obs ID 0200240101) 
on 2004 June 10, respectively. The European Photon 
Imaging Camera (EPIC), which consists of two MOS 
and one PN CCD arrays, was used with a medium fil- 
ter. The two EPIC/MOS cameras were operated in the 
Full-Frame Mode for a total exposure time of 16.5 ks 
for NGC 2346 and 19.2 ks for NGC 7026. The EPIG/pn 
camera was operated in the Extended Full-Frame Mode 
for a total exposure time of 11.1 ks for NGC 2346 and 
14.1 ks for NGC 7026. 

The XMM-Newton pipeline products have been pro- 
cessed using the XMM-Newton Science Analysis Software 
(SAS version 6.1.0) and the calibration files from the Cal- 
ibration Access Layer available on 2005 April 22. The 
event files have been screened to eliminate bad events 
and periods of high background. For the EPIC/MOS 
observations, only events with CCD patterns 0-12 (simi- 
lar to ASCA grades 0-4) were selected; for the EPIC/pn 
observation, only events with CCD pattern (single pixel 
events) were selected. To assess the background rate, we 
obtained light curves for each instrument in the 10-12 
keV energy range, where the counts are dominated by the 
background. The time intervals with high background, 
i.e., count rates > 0.25 cuts s^^ for the EPIC/MOS or 
> 1.0 cuts for the EPIC/pn, were discarded. The re- 
sulting exposure times for the EPIC/MOS and EPIG/pn 
observations were 16.5 ks and 10.9 ks for NGC 2346 and 
15.5 ks and 13.3 ks for NGC 7026, respectively 

3. RESULTS 

X-ray emission was clearly detected from NGC 7026 
but not from NGC 2346. We therefore describe the re- 
sults from NGC 7026 first and then determine an upper 
limit for the X-ray luminosity of NGC 2346. 

3.1. NGC 7026: Spatial Properties 

To measure the faint X-ray emission from NGC 7026, 
we used a 60" x 40" elliptical source aperture that en- 
compasses the entire optical nebula and a background 
aperture ~5 times as large as the source aperture lo- 
cated at a similar Y-coordinate of the detector. The 
EPIC/pn observation detected a total of 110 ± 15 
background-subtracted counts, corresponding to a count 
rate of 8.3x10-^ cuts s'^. The EPIG/MOSl and 
EPIC/M0S2 observations detected 26 ± 8 and 31 ± 8 
background-subtracted counts, corresponding to count 
rates of 1.7x10"'^ cuts s~^ and 2.0x10^^ cuts s~^, re- 
spectively. 

We have merged the EPIC/pn and EPIC/MOS images 
to increase the signal-to-noise ratio. The merged EPIC 
image is displayed in Figure la and a smoothed image 
in Figure lb. Diffuse X-ray emission is clearly detected. 
In order to compare the distribution of X-ray emission 
with the optical nebula, accurate alignment between X- 
ray and optical images is essential. We have determined 



the astrometric solution for the EPIC image with an ac- 
curacy of ~0'.'5, based on stellar sources that are detected 
in both this X-ray image and the Digitized Sky Survey. 
The relative alignment of the EPIC image and a Nordic 
Optical Telescope (NOT) [N II] image in Figure Ic is 
accurate to ~0'.'6. 

The distribution of diffuse X-ray emission from 
NGC 7026 is elongated in the North-South direction, 
matching the polar axis of the nebula. The position and 
orientation of the 50% contour of the smoothed EPIC 
image follow the bipolar lobes shown in the NOT [N ii] 
image (Fig. Ic). Furthermore, the X-ray emission shows 
two peaks that are projected within the bipolar lobes 
north and south of the central waist of NGC 7026. These 
excellent correspondences suggest that the diffuse X-ray 
emission from NGC 7026 is all confined within its bipo- 
lar lobes, similar to w hat is seen in the bipolar PN Mz 3 
()Kastner et al.ll2003(l . 

We have simulated the surface brightness of diffuse X- 
ray emission assuming the hot gas is distributed in two 
spherical shells with fractional shell thickness (AR/R) 
of 0.1, 0.2, 0.5, and 1.0. We find that the small number 
of counts and limited angular resolution precludes even 
a rough estimate of the geometry. Similarly, we cannot 
determine whether the local minimum along the waist is 
caused by higher circumstellar absorption. 

3.2. NGC 7026: Spectral Properties 

X-ray spectra of NGC 7026 have been extracted from 
the event files of the three XMM-Newton EPIC cam- 
eras using the same source and background apertures 
described in §3.1. Due to the limited number of counts 
detected by the EPIG/MOS observations, we will limit 
the spectral analysis to the EPIG/pn spectrum shown in 
Figure 2. This spectrum is soft, peaking at ~ 0.5 keV 
with a nearly flat distribution from 0.2 to 0.4 keV. Emis- 
sion between 0.7 and 1.1 keV is detected but with low 
S/N ratio. There is no detectable emission at energies 
greater than 1.25 keV. 

The EPIG/pn spectrum of NGC 7026, similar to other 
PNe, is suggestive of thermal plasma emission with its 
broad peak at ^0.5 keV corresponding to emission from 
the He-like triplets of N vi at ~0.43 keV and O vii at 
~ 0.57 keV. In order to model this X-ray spectrum, we 
have adopted the MEKAL thin plasma emission model 
([Kaa stra & Mewc 1993; Liedahl, Ostcrhcld, & GoldsteiiJ 
119951) . Spectral fits were carried out by folding the thin 
plasma emission model spectrum through the EPIG/pn 
response matrix, and comparing the modeled spectrum 
to the observed EPIC/pn spectrum in the 0.25-1.5 
keV energy range after the spectrum was resampled 
to energy bins with at least 16 counts. The best-fit 
was judged based on the statistics. It was found 
that the spectral fits did not accurately constrain the 
absorption column density, Nn, if it was allowed to 
vary as a free parameter. We therefore held A'h fixed 
at 3.2x10^^ cm^^, estimated from the observed loga- 
rithmic extinction in the H/3 line c(H/3) = 0.8 ± 0.1 
l|Kwitter fc Henrvll2001l: iHvung fc Feibelmanll200l and 
the canonical gas-to-dust ratio of Nyt/Ei B — V) = 
5.8 x 10^"'^ atoms cm^^ mag^^ (Bohlin. Savag ed Drake! 
Il978j) . The X-ray-emitting gas consists of shocked fast 
wind and mixed- in nebular gas. As the wind abun- 
dances determined from UV spectroscopic observations 
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Fig. 1.— XMM-Newton and NOT [N n] images of NGC 7026. Panel (a) displays the XMM-Newton EPIC raw image in the 0.4-1.1 keV 
band with pixel size 2'.'0. This energy band includes most of the photons detected from NGC 7026 as illustrated by its EPIC/pn X-ray 
spectrum shown in Fig. 2. Panel (b) shows a grey-scale presentation of an adaptively smoothed EPIC image using Gaussian profiles with 
FWHM ranging from I'.'S to 4". The X-ray contours overplotted correspond to 50%, 75%, and 95% of the peak level. Panel (c) displays 
the NOT [N ll] image overplotted by these same X-ray contours. Panel (d) shows a color-composite image comprised of the smoothed 
XMM-Newton EPIC image (blue) and HST images obtained with the F502N (green) and F658N (red) filters to show [O III] and [N II] 
emission. 



l|Koesterkdl2nnit IHa mann. Todt fc Graefeneill2?ffll dif- 
fer from the nebu l ar abundances determined from optical 
spectra ()K witter fc Henrvll2001h . we have made spectral 
fits using wind or nebular abundances separately. 

The first spectral fits use the nebular abundances, 
with the number ratio of He, O, N, Ne, Ar, CI, 



and S to hydrogen being 0.14, T.SxlO"*, S.SxlO"", 
1.8x1 0-^, 3.1x10-6, 3.7x10- ^, and 1.4xl0-^ respec- 
tively l|K witter fc Henrvll200lD . The corresponding neb- 
ular abund ances of these eleme nts relative to the so- 
lar values ijAnders fc Grevessd 'l989) are 1.43, 0.86, 
4.9, 1.8, 1.5, 2.0, and 0.96, respectively. Solar abun- 
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Fig. 2.— (left) XMM-Newton EPIC/pn background-subtracted spectrum of NGC 7026 overplotted with the best-fit MEKAL model 
assuming nebular abundances (solid line), (right) the same spectrum overplotted with the best-fit MEKAL model assuming the abundances 
of the central star. For plotting purposes, both the spectrum and model are binned in 50 eV wide energy bins to preserve the EPIC/pn 
spectral resolution of ~ 100 eV at 1 keV. Beneath each spectrum the residuals from the best-fit are plotted. 



dances are used for the other elements, which have low 
abundances and are not expected to contribute signifi- 
cantly to the diffuse X-ray emission. We have also as- 
sumed solar abundances for the foreground interstellar 
abs orption, and adopted absorptio n cross-sections from 
iBahicinska-Church fc McCammonI l(1992il . The result- 
ing best-fit model with nebular abundances and a fixed 
A^H is overplotted on the EPIC/pn spectrum in Fig- 
ure 2a. The plasma temperature of the best-fit model 
is T (l.l±oi) X lO^ K (i.e., kT = 0.095toml keV). 

The second spectral fits use stellar abundances. The 
central star of NGC 7026 has been classified as a 
hydrogen-deficient WC star with the mass-fractions of 
He, C , and O being 70, 20, and 10%, respectively 
ijKoest erkc 2001; Hama nn et ahfe OOSl. In order to make 
spectral fits using a MEKAL model we have assumed a 
He/H ratio of 200 and scaled the abundances of C and 
O according to their mass fractions. Thus, the values 
adopted for the stellar wind abundances relative to the 
solar values were 20, 500, and 80 for He, C, and O, re- 
spectively. All other abundances were set to zero and the 
absorption column density, A^h, was again held fixed with 
a value of 3.2 x 10^^ cm~^. The resulting best-fit model is 
shown in Figure 2b. This model with stellar abundances 
has roughly the same temperature as the best-fit using 
nebular abundances and a slightly higher reduced x^- 
The most obvious difference between the best-fit models 
using nebular and stellar abundances is that the model 
with stellar abundances underestimates the emission at 
~0.5 keV due to the absence of N in the stellar chemi- 
cal composition. Clearly the quality of both spectral fits 
are limited by the relatively small number of counts in 
the X-ray spectrum and it is not possible to convincingly 
distinguish whether the chemical composition of the X- 
ray emitting gas in NGC 7026 better matches that of the 
central star or the nebula. 

From the EPIC/pn spectrum we can also deter- 
mine that the volume emission measure is 3.6lQ g x 
1055(__d_)2 ^^-3^ ^Yiere d is the distance to NGC 7026 
in pc. The distance has been estimated to be 



1700;^^ii pc based on th e extinction to NGC 7026 
ijHvung fc FeibelmanI l2004fl and o n the e xtinction- 
distan ce relationship derived by iSolf fc W cinbcrgeil 
ifTosl . The observed (absorbed) X-ray flux of NGC 7026 
in the 0.2-2.5 keV energy range is (8.8 ± 1.2) x 10"^^ 
ergs cm^^ s^^, and the intrinsic (unabsorbed) X-ray 
flux is {l-ito'D x 10-12 grgs cm-2 s'^. The total X- 
ray luminosity in the 0.2-2.5 keV energy range is Lx = 
(4.5l?:°) X 1032(^)2 ergs s'^. 

3.3. NGC2346 

The XMM-Newton observations of NGC 2346 are com- 
pared with a Hubble Space Telescope {HST) WFPC2 
[N ii] image in Figure 3. There is no statistically sig- 
nificant X-ray emission associated with either the nebula 
or the central star in these observations. The EPIC/pn 
observations have been used to derive an upper limit 
for its X-ray luminosity using a source aperture encom- 
passing the optical nebula. The 3a upper limit for the 
EPIC/pn count rate in the 0.3-1.25 keV band is 4.5 x 10"^ 
cuts s^i. To derive an upper limit of the intrinsic (un- 
absorbed) flux we assume a thin plasma emission model 
with similar temperature and nebular abundances to the 
one that describes the X-ray emission from NGC 7026 
(see §3.2). The foreground absorption is estimated based 
on the obser ved logarithmic extinction at the H/3 line, 
c(H/3)=0.75 ijCiarduUo et al.lll999j) . which corresponds 
to a column density of Nyf= 3.QxlO^^ cm^ '^. At a dis- 
tance to NGC 2346 of 690 pc l)Terzianlll99^ . the 3 -a up- 
per limit of the unabsorbed flux corresponds to an X-ray 
luminosity <4xl0'^i ergs s~^. 

4. DISCUSSION 

The XMM-Newton observations of NGC 7026 clearly 
show diffuse X-ray emission while the observations of 
NGC 2346 do not detect significant X-ray emission. We 
will now use these observations to discuss the production 
and evolution of hot gas in bipolar PNe. 

4.1. NGC 7026 versus NGC 2346 
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Fig. 3.— XMM-Newton and i?5T [N ll] images of NGC 2346. (left) XMM-Newton EPIC image in the 0.25 - 1.5 keV band with pixel 
size 2''0. The EPIC/pn and EPIC/MOS images have been merged and the resulting image smoothed using Gaussian profiles with FWHM 
4" to increase the signal-to-noise ratio. The X-ray contours overplotted corresponds to the 1, 3, 5, and 9-(T over the background, {right) 
HST WFPC2 [N ll] image overplotted by the same X-ray contours. 



The spatial distribution of the X-ray emission in 
NGC 7026 suggests that the X-ray-emitting gas is con- 
fined within the bipolar lobes, but its spatial distribution 
is unknown. Therefore, we define the hot gas filling fac- 
tor, e, as the ratio of the emitting volume to the total neb- 
ular volume and use a reference value of 0.5 in our analy- 
sis. Using the volume emission measure determined from 
the spectral fits (in §3.2) and approximating the emitting 
volume by a 40" x 15" x 15" rectangular box, the rms elec- 
tron density is ~ 22 (d/1700pc)'i/2(£/o.5)-i/2 ^m-^. 
This density and the best-fit plasma temperature can in 
turn be used to derive the thermal pressure of the hot gas, 
Pth ~ 1.9 X UckT ~ 6 X 10-9 (d/1700pc)-i/2 
dynes cm"^, assuming H and He are fully ionized and 
that the number ratio of He/H is 0.14 (the nebular abun- 
dance) . 

We can compare this hot gas pressure to the thermal 
pressure of the warm nebular shell of NGC 7026 to see 
whether the nebular expansion is driven by th e ther- 
mal pressure as expected by the iWeaver et all (^23) 
mo del. The electron density and temperature reported 
by iKwitter fc H enrv^ ("2001") are generally lower than 
those reported by Hyung & Fcibclman ( 2004) . The low- 
est temperature and density are derived fro m the [S II ] 
lines, 6,400 K and 3,300 cm'^ (jKwitter fc Henrvll200H) . 
resulting in a thermal pressure of 6.3 x IQ-^ dynes cm"^. 
The temperatures and densities derived from a variety 
of fines by Hvung & Feibelman (2003) are ~ 9,000 K 
and ~ lO'* cm"'^, resulting in a thermal pressure of 
~ 2.3 X IQ-® dynes cm"^. These measurements were 
made at the bright nebular waist, where the density is 
expected to be the highest. While the thermal pressure 
of the nebular gas is higher than that of the hot gas in the 
waist region, it is likely that the pressure of the nebular 
gas in the lobes is comparable to that of the hot gas. 

NGC 2346 is a larger and more evolved bipolar PN 
than NGC 7026. No fast wind from the central star of 
NGC 2346 has been detected, as the lUE observations 
of the stellar C IV and N V lines do not show P Cygni 



profiles l|Patriarchi fc Perinottoll991|) . If the central star 
of NGC 2346 had a fast stellar wind in the past, the hot 
gas resulting from the shocked fast wind has either cooled 
or is too rarefied to be detected. 

While the presence of a strong fast stellar wind is a 
prerequisite for the production of detectable diffuse X- 
ray emission, the lobe structure of bipolar PNe may 
also play an essential role. NGC 2346 has an open-lobe 
structure and lacks detectable diffuse X-ray emission and 
NGC 7026 has a closed-lobe structure and detectable X- 
ray emission. Previous X-ray observations of bipolar PNe 
have detected only two other nebulae, NGC 7027 and 
Mz3, and both have closed-lobe structures. It is possible 
that open lobes allow shocked- wind to be distributed to a 
larger volume and the lower density precludes the detec- 
tion of X-ray emission. In the case of NGC 2346, hot gas 
is not produced currently because of a lack of fast stel- 
lar wind. Even if it had a fast wind in the recent past, 
its open lobes were not able to retain the shocked stel- 
lar wind to prolong the detectability of its diffuse X-ray 
emission. 

X-ray emission is not detected from the central star 
of either NGC 2346 or NGC 7026. The central star of 
NGC 2346 is a hot white dwarf with an A5 companion. 
A hot white dwarf's photospheric emission can reach the 
soft X-ray energy range, but is easily absorbed by the 
nebular gas or foreground interstellar material. A-type 
stars have neither fast stellar winds nor active coronae, 
and therefore are not expected to be X-ray sources. Even 
if the white dwarf possesses an accretion disk, the tem- 
perature of the accretion disk can not be high enough to 
produce X-ray emission. NGC 7026 has a WC-type cen- 
tral star. No X-ray emission has ever been detected from 
a single WC star even among the massive Wolf-Rayet 
stars (Oskinova et al. 2003) because of the high opacity 
in their stellar atmospheres. 
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To date, only three bipolar PNe have been detected 
m X-rays. Among these, Mz3 has a symbiotic central 
star and its X-ray emission appears to originate from 
interac tions between its bip olar jets and the nebular ma- 
terial (KastMriEiOllDOS). The other two, NGC 7026 
and NGC 7027, are bona fide bipolar PNe. The X-ray 
luminosity of NGC 7026 is one of the highest among 
all PNe detected in X-rays. Its hot gas mass, '^0.002 
(j7^^)^/^(e/O.5)^/^M0, is also high compared to other 
PNe. Note however that these values may be over- 
estimated because the X-ray-emitting lobes may have 
lower circum-nebular absorption than the waist region 
where the extinction was measured. The young bipolar 
PN NGC 7027 has a lower X-ray luminosity and smaller 
hot gas mass, but a higher hot gas temperature and 
density ijKastner et alJl2001|) . The differences between 
NGC 7026 and NGC 7027 might be due to an evolution- 
ary effect, but more positive X-ray detections of bipolar 
PNe of intermediate ages are needed to investigate the 



properties and evolution of their hot gas. 

Chandra has observed five other bipolar PNe but none 
were detected. Based on these non-detections we sug- 
gest that future X-ray observations of bipolar PNe should 
avoid objects that have open bipolar lobes and objects 
with low velocity outflows (<300 km s~^). Instead, ob- 
servations of objects with closed bipolar lobes, to confine 
hot gas, and with kinematic ages of no more than a few 
thousand years would appear best suited for understand- 
ing the generation and evolution of hot gas in bipolar 
PNe. 
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